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OBJECT 


Studies  toward  the  establishment  of  a theory  of  recoilless  rifle  blast. 


SUMMARY 

The  study  of  the  blast  from  recoilless  rifles  requires  a theoretical  analysis  as 
a preliminary  step  to  an  extensive  experimental  mapping  of  the  blast  field.  To  this 
end,  a theoretical  analysis  of  the  flow  from  a nozzle  has  been  made.  A single,  wedge 
shaped  nozzle  with  7°  half-angle  of  divergence  was  considered.  The  ratio  of  exit 
height  to  entrance  height  of  the  nozzle  is  V2: 1 . Chamber  pressure  and  temperature 
have  been  assumed  to  be  5400  psi  and  2400°  K,  respectively.  The  principal  feature  of 
the  flow  is  a concentrated,  almost  stationary,  region  of  strong  shock  called  the  Mach 
disk.  This  was  not  located  in  the  course  of  the  usual  calculations  because  of  the 
limited  assumptions  of  the  theory.  On  the  basis  of  the  growth  of  turbulence  in  the 
boundary  layer,  the  distance  of  the  Mach  disk  from  the  exit  section  of  the  nozzle  was 
found  to  be  between  40  and  140  exit  diameters,  the  first  value  being  the  more  probable. 
The  pressure  and  temperature  after  the  shock  are  40  psi  and  2300°  K in  the  first  case, 
14.6  psi  and  3100°  K in  the  second  case.  Mach  numbers  as  high  as  ten  occur  in  the 
flow.  A detailed  map  of  the  flow  up  to  a point  140  exit  diameters  away  from  the  exit 
section  is  given.  This  map  should  provide  a guide  for  the  study  of  the  blast  from  re- 
coilless rifles,  since  it  enables  cne  to  interpret  better  the  results  of  experimental 
studies.  Furthermore,  the  location  of  the  Mach  disk  can  serve  as  a criterion  for 
further  improvements  in  the  theory. 
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INTRODUCTION 


The  study  of  the  blast  from  recoilless  rifles  must  proceed  in  two  parallel  modes, 
the  experimental  and  the  theoretical.  The  first  step  in  the  theoretical  analysis  of 
the  blast  of  gas  from  recoiJless  rifles  is  the  study  of  flow  from  a very  high  pressure 
nozzle  fchamlter  pressure  greater  than  1000  psi).  This  report  gives  the  first  approxi- 
mation in  the  study  of  the  free  gas  jet  coming  from  such  a nozzle.  Such  a study  is 
only  the  first  of  a series  of  calculations,  but  by  itself  it  yields  significant  infor- 
mation about  high  pressure  jets  in  general,  and  thus  may  be  useful  an  a guide  in  long 
range  planning. 

It  is  a general  feature  of  high  pressure  jets  that  they  are  separated  into  two 
distinct  regions  by  a very  strong  shock  pe rpenil  i cu 1 a r to  the  axis  of  the  nozzle, 
called  the  principal  Mach  di.-k.  Studies  of  the  region  downstream  fi.e.,  the  region 
farthest  from  the  nozzle)  of  this  shock  may  be  done  independently  of  work  on  the  up- 
stream side. 

This  report  ''vers  the  region  of  flow  fin  the  two-dimensional  case)  lie  een  the 
nozzle  and  th.  principal  Mach  disk.  The  chief  results  of  this  wor1:  are  contained  in 
the  maps  of  the  flow  field  given  in  Figures  1 and  2 and  Table  I.  These  give  approxi- 
mate values  of  pressure  and  velocity  at  any  point  in  the  indicated  region  of  gas  flow. 
Preparations  are  being  made  to  correlate  the  theoretical  analysis  of  the  blast  with  an 
experimental  map  of  the  flow  field. 


BASIC  DEFINITIONS 


The  specialized  vocabulary  of  fluid  dynamics  requires  a set  of  definitions  for 
the  full  under standing  of  any  report  on  some  phase  of  that  subject. 


Sonic  Velocity 

The  velocity  of  sound  in  the  medium  at  the  point  under  consideration;  usually  as- 
sumed to  be  the  square  root  of  the  ratio  of  the  specific  heats  of  the  gas  times  the 
local  pressure  divided  by  the  local  density. 


Hodograph 

Refers  to  the  "velocity"  plane,  the  plane  whose  coordinates  are  the  components  of 
the  gas  velocity  vector.  Bv  the  hodograph  ’transformation,1  the  determinat ion  of  the 
flow  i.t  a point  in  the  physical  (x,y)  plane  is  made  by  examining  the  flow  at  a corre- 
sponding point  oi  the  hodograph  plane. 
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Sonic  Circle 


to 


The  locus  of  those  points  in 
the  local  sound  velocity.  This 


the  hodograph  plane  where 
locus  is  a circle. 


the  gas  velocity 


is  eoual 


H.ach  Number 

A dimensionless  parameter  equal  to  the  ratio  of  gas  velocity  to  local  sound 
velocity.  In  supersonic  flow  the  Mach  number  is  al-.vays  greater  than  one.  On  the 
sonic  circle  the  Mach  number  is  one. 


Hach  Angle 

The  angle  whose  sine  is  the  reciprocal  of  the  Mach  number.  This  angle  is  real 
for  supersonic  flows  only. 

Character i sties 

A network  of  lines  which  is  derived  from  the  equations  cf  motion  and  is  used  to 
characterize  a particular  solution.  The  analytic  definition  is  given  in  Appendix  A. 
In  two  dimensions  there  are  two  characteristics  through  each  point.  The  flow  Quanti- 
ties (pressure,  density,  velocity)  arc  determined  at  any  point  by  their  values  on 
tho,.e  neighboring  points  which  are  on  one  of  the  two  characteristics  passing  through 
the  point  in  question,  if  the  characteristics  air  real.  This  is  the  case  in  super- 
sonic flow. 


Rotat  i octal  ity  and  Vortici  ty 

These  quantities  refer  to  the  angular  momentum  present  locally  in  a flow  without 
singularities.  The  vorticity  is  proportional  to  the  local  angular  momentum;  when  the 
vorticity  is  zero  in  a region,  the  flow  throughout  the  region  is  said  to  be  lrrota- 
t lonal . 

Shock  Bottle 

A configuration  of  shocks  usually  seen  in  high  pressure  jets,  such  as  are  found 
at  the  muzzles  of  guns.  There  are  two  shocks  near  the  jet  fx>undary,  curving  inward  as 
one  goes  away  from  the  exit  of  the  nozzle,  terminating  in  a Mach  disk.  This  structure 
has  the  outlines  of  a bottle,  hence  the  name.  There  are  usually  two  snocks  leaving 
the  Mach  disk  which  look  like  reflections  of  the  two  shocks  forming  the  sides  of  the 
bottle.  The  two  shocks  nearest  the  nozzle  exit  are  called  intercepting  shocks.  Those 
on  the  downstream  side  of  the  Mach  disk  are  called  reflected  shocks.  (Figure  4). 
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Boundary  Layer 


A very  thin  region  of  transition  in  what  is  usually  considered  the  boundary  of 
the  flow.  The  flow  in  this  region  is  considered,  in  general,  as  being  determined  by 
the  flow  in  the  main  portion  of  tne  stream,  i.e.,  the  flow  in  which  no  boundary  lay0*- 
is  taken  into  account.  Tne  boundary  layer  is  subsonic  and  often  turbulent. 


BASIC  ASSUMPTIONS 

The  following  assumptions  are  made: 

The  most  general  feature*-  of  this  type  of  flow  are  contained  in  the  flow  from  a 
circular  nozzle  whose  expansion  ratio  and  divergence  are  small.  This  assumption  is 
not  questionable,  except  when  one  wishes  to  consider  the  flow  on  the  axis  near  the 
nozzle.  High  pressure,  in  general,  tends  to  decrease  the  dependence  of  the  flow  on 
nozzle  shape  at  l<irge  distances  and  increase  it  at  small  distances  from  the  exit 
section.  The  critical  distance  is  of  the  order  1 L,  where  L is  the  principal 
nozzle  dimension  and  M is  the  Mach  numlier  at  the  exit  section. 

The  flow  is  isentropic  from  the  throat  of  the  nozzle  to  the  first  shock  reach  d 
in  t.ie  atmosphere.  This  neglects  the  effects  of  minor  shocks  arising  in  the  nozzle 
from  surface  defects,  and  shocks  on  the  axis  from  lip  effects,  both  of  which  must  be 
neglected  if  the  calculations  are  to  be  made  in  a reasonable  time.  Conversations  with 
various  authorities  in  the  field  reveal  that  the  first  of  these  effects  is  not  likely 
to  be  large,  while  the  last  can  not  be  estimated  at  this  time. 

The  flow  is  i rrotat ional ; the  ratio  of  specific  heats  is  constant;  secondary  com- 
bustion can  be  neglected.  These  assumptions  are  necessary  to  facilitate  computation. 
I r rot  a t ional i ty  is  associated  with  isentropic  flow.  Variations  in  the  ratio  of  spe- 
cific heats  can  be  neglected,  if  previous  work  in  the  field  is  any  indication.  The 
neglect  of  secondary  combustion  is  necessary  in  order  to  compute  the  points  where 
secondary  combustion  can  occur.  In  other  words,  such  combustion  is  viewed  as  a small 
perturbation  of  the  flow,  like  the  boundary  layer. 

Boundary  layer  effects  are  to  be  excluded.  This  assumption  is  usually  made  in 
calculations  like  these,  and  indeed  it  is  necessary  to  do  so;  but  this  question  will 
be  treated  more  fully  in  the  body  of  the  report. 

Sufficient  accuracy  is  obtained  when  the  flow  is  determined  via  the  method  of 
characteristics  for  the  corresponding  two-dimensional  case.  This  is  standard  pro- 
cedure in  computing  axially  symmetric  flows.  There  are  points  where  the  correction  is 
likely  to  be  considerable;  but  in  view  of  the  various  other  approximations  it  is 
doubtful  that  the  computations  necessary  to  make  this  correction  are  justified  at  this 
time.  The  two-dimensional  approximation  will  be  treated  more  fully  in  the  body  of  the 
report . 
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PROCEDURE 


A two-dimensional,  wedge  shaped  nozzle,  whose  expansion  ratio  is  ^2: 1 , and  whose 
half-angle  of  divergence  is  7°,  is  chosen.  The  chamber  pressure  is  taken  as  5400  psi. 
The  chamber  temperature  is  taken  as  2400°  K.  Atmospheric  pressure  is  taken  as  14.7 
psi,  y = ratio  of  specific  heats  = 1.4. 

The  method  of  characteristics  is  used  to  map  the  corresponding  two-dimensional 
flow  field.  Tile  characteristics  method  is  adequately  described  in  standard  texts*  and 
is  reported  more  fully  in  Appendix  A.  Oily  a general  description  of  the  method  will 
lie  given  here. 

The  flow  equations,  which  are  discussed  in  Appendix  A,  are  partial  differential 
equations  of  the  hyperbolic  type  (in  the  case  of  supersonic  flow).  The  construction 
of  the  solution  of  a partial  differential  equation  of  hyperliolic  type  was  formulated 
in  great  generality  by  Kiemann  around  1860.  Earlier  the  method  had  been  suggested  for 
electrodynamics  by  Huyghens , on  physical  grounds.  Thus,  it  is  often  the  practice  to 
refer  to  the  Riemnrn  method  of  solution  as  a Huyghens  construction,  even  though  the 
equation  is  not  from  electrodynamics. 

In  general,  linear  equations  permit  of  much  more  el egan t and  soph i s t i ca ted  metliods 
of  solution,  and  the  method  of  characteristics  is  reserved  for  such  tasks  as  de- 
limiting the  region  in  which  a solution  exists.  This  is  particularly  the  case  when 
the  coefficients  are  constant;  hut  if  the  coefficients  vary  widely,  and/or  the  equation 
is  nonlinear,  such  methods  fail  and  recourse  must  be  had  to  the  characteristics  method. 

The  success  of  this  method  is  limited  to  two  dimensions,  by  and  large,  by  compu- 
tational difficulties.  The  accuracy  is  limited  by  the  size  of  the  mesh,  and  even 
though  the  initial  steps  are  quite  small,  the  mesh  becomes  very  large  if  there  is 
larp-'  expansion,  and  there  is  a resulting  loss  in  accuracy.  One  must  assume  the  state 
of  flow  along  some  initial  curve.  Boundary  conditions,  such  as  containment  within 
walls  and  the  existence  of  a free  surface  in  the  atmosphere,  can  be  applied  auto- 
matically. Other  types  of  boundary  conditions  ure  difficult  to  apply  but  luckily  do 
not  often  arise. 

It  should  be  understood  that  the  method  is  valid  only  up  to  the  formation  of  a 
shock,  but  it  can  be  continued  formally  beyond  the  shock  to  give  a lower  order  of 
approximation  there. 

Appendix  A describes  the  analytical  techniques  for  obtaining  pressure  and  density 
from  a map  of  the  flow  characteristics.  The  details  of  preparing  such  a map  are  given 
in  Appendix  B. 


*A.  Ferri,  "Elements  of  Aerodynamics  of  Supersonic  Flows,"  Macmillan,  New  York,  1949; 
pp  19-40  (ch  2). 
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RESULTS  AND  DISCUSSION 


i 


Figure  1 is  the  map  of  the  flow  in  the  nozzi  , and  of  the  flow  immediately  beyond 
the  nozzle  exit.  Figure  2 is  the  map  of  the  entire  flow  up  to  the  principal  Mach 
disk.  Figure  3 is  a plot  of  the  static  pressor^  along  the  line  of  fire,  behind  the 
gun.  Tables  I and  II  enable  one  to  find  the  pressure,  density,  velocity,  etc.,  in  any 
particular  square  of  the  flow  map. 


Flow  in  the  Nozzle 

Flow  in  the  nozzle  is  almost  identical  with  that  predicted  by  hydraulic  theory. 
This  will  not  be  the  case  in  general  for  nozzles  of  much  higher  angles  of  divergence. 
The  most  striking  feature  of  the  flow  immediately  beyond  the  exit  section  is  the  sharp 
turn,  through  57°  (Figure  1),  made  by  the  boundary  streamline  in  order  to  form  a free 
surface  at  a pressure  of  14.7  psi.  From  this  point  to  the  region  of  influence  of  the 
principal  Mach  disk  the  boundary  curves  back  towards  the  axis. 

Formation  of  the  Shock  Bottle 

Rarefaction  waves  from  the  exit  section  are  reflected  from  the  boundary,  anti  re- 
turn to  the  main  st  earn  as  compression  waves  which  form  an  envelope  called  the  inter- 
cepting shock  (the  sides  of  the  shock  bottle).  Computation  of  the  flow  is  continued 
formally  in  the  region  between  this  shock  and  the  boundary.  Neglect  of  shock  condi- 
tions gives  a lower  order  of  accuracy  there.  This  shock  is  comparatively  weak,  and 
corrections  are  not  large,  except  in  the  neighborhood  of  the  Mach  disk  and  its  re- 
flected shock. 

The  inclination  of  the  shock  to  the  stream  direction  is  almost  constant  at  24°. 
The  rot  a t iona 1 i ty  introduced  in  the  downstream  flow  depends  on  the  entropy  gradient 
there.  The  entropy  gradient  will  be  high,  since  the  variation  in  Mach  numbers  up- 
stream of  the  shock  is  from  6 to  10,  roughly. 


Pressures  in  Front  of  the  Shock  Bottle 

Static  pressure  along  the  axis  is  easily  plotted  and  the  results  are  given  in 
Figure  3.  Units  of  length  are  arbitrary  and  are  determined  by  the  half  height  of  the 
throat  which  is  1.25  units  of  length.  The  point  of  inflection  corresponds  to  a focus 
in  the  actual  three-diir.ensional  case,  where  the  sudden  change  in  the  streamline  slope 
at  the  exit  section  makes  itself  felt  as  a logarithmic  singularity  at  th'*  indicated 
point  on  the  axis.  That  point  is  the  intersection  of  cone  whose  base  is  the  exit 
section  and  whose  apex  angle  is  the  exit  Mach  angle.  This  effect  was  discovered  by 
S.  A.  Pai  at  the  Institute  for  Fluid  Dynamics,  University  of  Maryland. 
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Pressure  along  other  radial  lines  requires  extensive  interpolation  before  it  can 
be  tatxilated.  The  very  low  pressures  (and,  consequently,  temperatures)  and  very  high 
Mach  numbers  produced  in  the  flow  should  be  noted.  These  point  to  the  possibility  of 
condensation  and  consequent  condensation  shocks. 

Location  of  the  Mach  Disk 

It  was  hoped  that  the  method  of  calculation  employed  here  would  yield  the  loca- 
tion and  extent  of  the  principal  Mach  disk.  It  was  surprising  to  find  that  this  did 
not  occur.  The  intercepting  shock  reflected  from  the  symmetry  axis  and  it  was  not 
possible  to  discern  any  reason  for  the  location  of  a Mach  disk.  The  shock  pattern 
found  was  revealed  to  f>e  well  within  critical  limits  and  on  the  basis  of  this  work 
alone  there  was  no  1 eason  to  expect  a Mach  disk. 

Nevertheless,  such  disks  have  invariably  been  found  in  jets  whose  pressures, 
while  more  than  a few  atmospheres,  were  much  lower  than  those  used  here.  Moreover, 
the  shock  bottle,  whose  top  is  a Mach  disk,  has  been  observed  time  and  again  at  the 
muzzles  of  guns.  D.  C.  Pack  of  the  1:  stitute  for  Fluid  Dynamics,  University  of 
Maryland,  has  said  that  similar  difficulties  had  been  noted  by  a number  of  his  stu- 
dents. Conseauent ly , it  is  believed  that  some  other  mechanism,  possibly  one  hereto- 

ry 

fore  purposely  neglected,  is  responsible  for  this  phenomenon.  Courant  and  Friedrichs” 
express  the  opinion  that  the  process  of  mixing  with  the  surrounding  air  is  largely  the 
determining  factor. 

Some  hypotheses  about  the  formation  of  the  Mach  disk  are: 

1.  The  disk  is  a stable  conf igui at  ion  established  early  in  the  formation  of  the 
jet  by  the  action  of  transient  phenomena. 

2.  In  steady  state  flows,  the  transition  of  the  boundary  layer  from  laminar  to 
turbulent  flow  causes  a disturbance  in  the  main  stream  which  sets  up  a Mach  disk. 

In  either  case  the  mechanism  is  the  same:  vorticity  is  produced  at  some  point  or 
region  and  propagates  into  the  main  stream  (via  the  heat - conduc * ion- convect ion  equa- 
tion). Extra  vorticity  can  be  introduced  into  the  flow  only  by  curving  the  shocks 
more  strongly.  After  a point  an  oblique  shock  is  no  longer  possible  and  a normal 
shock  (Mach  disk)  comes  into  being.  Presumably,  thrre  is  a point  where  no  laminar 
flow  is  possible  anywhere  in  the  jet  and  the  usual  jet  structure  breaks  down. 

In  the  first  case,  vorticity  is  produced  either  by  mixing  with  the  air  or  inter- 
ference from  the  nozzle  lips.  Then  at  some  point  in  the  early  development  of  the  jet, 
the  vortex  ring  and  the  back  pressure  associated  with  the  short  length  of  the  growing 


^Courant  and  Friedricks,  "Supersonic  Flow,"  Interscience,  Nf  » York,  p.  392. 
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j e t causes  the  flow  to  choo'"  one  possible  pattern  rather  than  another,  the  chosen 
oattern  being  stable.^  Such  ambiguities  in  supersonic  flow  are  well  known. 

In  the  second  case,  the  boundary  layer  changes  abruptly  from  laminar  to  turbulent 
flow.  In  such  a case  the  pressure  on  the  main  !>ody  of  the  stieam  changes  sharply  as 
the  pressure  gradient  in  the  boundary  layer  changes  from  that  characteristic  of 
laminar  flow  to  that  characteristic  of  turbulence.  Such  a disturbance  is  propagated 
along  a characteristic  line  to  the  rest  of  the  flow.  Such  a disturbance  can  act  not 
only  as  a source  of  vorticity  but  also  as  a strong  pressure  pulse.  When  this  charac- 
teristic hits  the  intercepting  shock,  the  curvature  of  the  shock  must  change  in  order 
to  take  this  new  disturbance  into  account.  If  the  pressure  is  sufficiently  high  (the 
critical  value  depends  upon  *ho  nozzle)  the  shock  must  he  normal.  Unfortunately,  the 
theory  of  the  boundary  layer  transition  is  not  formulated  sufficiently  at  present  for 
definite  predictions  to  be  made.  Recent  work  on  the  transition  from  laminar  to  turbu- 
lent boundaiy  layer  flows  shows  that  this  transition  is  associated  with  Reynolds’ 
numbers  of  the  order  U)'J  for  conical  flow  (using  velocity  along  the  surface),  106  for 
flow  over  a flat  plate.  A quick,  rough  calculation  shows  that  this  magnitude  Reynolds’ 
number  is  likely  to  be  reached  at  a point  one-third  of  the  way  along  the  boundary  of 
the  jet  between  the  exit  section  and  the  point  of  maximum  expansion.  This  will  reduce 
the  distance  of  the  Mach  disk  from  the  exit  section  to  about  one- third  of  the  value 
shown  in  Figure  2. 

It  is  reasonable  to  assume  that  the  boundary  layer  transition  will  occur,  if  at 
all,  between  the  lips  of  the  nozzle  and  the  point  where  t lie  boundary  becomes  hori- 
zontal. The  first  point  is  the  start  of  mixing  with  the  surrounding  atmosphere,  the 
second  the  only  other  point  distinguishable  on  physical  grounds  from  any  other  point 
on  the  boundary.  In  other  words,  the  assumption  is  made  that  the  transition  occurs  in 
the  first  half  of  the  main  jet  structure.  This  assumption  is  justified  by  Schlieren 
photograohs  of  jets  available  from  a number  of  sources/ 

For  the  purposes  of  this  analysis,  the  point  of  intersection  of  the  characteristic 
from  the  horizontal  part  of  the  boundary  with  the  intercepting  shock  was  chosen  as  the 
location  of  the  normal  shock.  This  represents  the  maximum  distance  away  from  the 
nozzle  of  the  Mach  disk.  This  distance  is  of  the  order  of  140  exit  diameters.  There 
are  many  indications  that  this  is  two  or  three  times  too  large,  indicating  a tran- 
sition point  somewhat  closer  to  the  nozzle.  Such  an  error  is  caused  not  only  by  im- 
proper location  of  the  critical  characteristic  line,  but  also  by  neglect  of  any  shocks 
or  similar  details  on  the  axis  caused  by  nozzle  details.  Otherwise,  errors  of  this 
magnitude  are  not  expected. 


*Har  tmann,  J and  Lazarus,  F. , "The  Air-jet  with  a Velocity  Exceeding  That  of  Sound." 
Phil.  Mag.  31,  7,  35-50,  1941. 

^Hartmann  and  Lazarus,  loc.  c i t . : Paymnn  and  Sheperd,  "Explosion  Waves  and  Shock 

Waves,"  Proc.  Roy.  Soc.  A.  186,  plate  17,  Figures  10-20,  (1946). 
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Experimental  studies  of  the  flow  of  powder  gas  through  a conical  nozzle  similar 
in  cross-section  to  the  one  studied  here  show  that  the  Mach  disk  occurs  about  15  exit 
diameters  from  the  nozzle  exit.  This  corresponds  to  the  point  (Figure  2)  where  the 
characteristic  leaving  from  the  focus  intersect?  the  boundary.  As  a numf>er  of  authori 
ties  have  mentioned,  this  corresponds  to  the  first  point  where  the  entire  flow  is 
aware  of  the  presence  of  the  atmosphere.  Hence,  this  represents  the  point  closest  to 
the  nozzle  that  the  Mach  disk  can  occur.  This  can  then  be  used  as  a criterion  for 
locating  the  disk  with  considerably  better  accuracy  than  the  method  mentioned  previ- 
ously. These  experimental  studies  were  conducted  some  time  after  the  main  Ixxly  of  the 
report  was  written. 


CONCLUSIONS 

The  study  presented  here  represents  an  attempt  to  understand  the  most  general 
features  of  the  jet  from  a nozzle  backed  by  very  high  pressure.  As  soon  as  a reliable 
method  for  locating  the  Mach  disk  is  available,  this  study  will  provide  an  adequate 
representation  of  the  flow  in  the  section  nearest  the  nozzle.  In  the  meantime  this 
study  establishes  the  maximum  extent  of  that  flow. 

The  pressures  indicated  in  the  graph  in  Figure  3 and  the  corresponding  dynamic 
pressures  give  a good  indication  of  the  forces  to  be  expected  in  the  region  upstream 
of  the  Mach  disk.  In  the  downstream  region  an  analysis  of  the  mixing  with  the  air  is 
needed  before  similar  conclusions  can  l>e  drawn. 

The  extrapolation  to  three  dimensions  can  be  carried  out  formally  by  a method 
outlined  in:  R.  Sauer,  "Method  of  Characteristics  for  Three-Dimensional  Axially  Sym- 
metric Super-  nic  Flows,"  N.A.C. A.  T.M.  1133,  Washington,  1947.  Examples  given  in 
this  reference  indicate  that  linear  dimensions  are  decreased  by  a few  per  cent,  in 
general.  This  arises  from  the  extra  room  for  expansion  provided  by  the  third  degree 
of  freedom.  More  quantitative  statements  than  this  can  not  l>e  made  at  this  time  be- 
cause of  the  difficulty  in  locating  the  principal  Mach  disk.  The  only  effect  opposing 
the  decrease  in  linear  dimensions  is  that  of  secondary  combustion.  rhe  low  tempera- 
tures and  pressures  before  the  disk  and  the  high  temperatures  and  pressures  after  the 
disk  indicate  that  most  secondary  combustion  will  occur  after  the  disk. 

It  should  be  pointed  out  that  even  a phenomenological  theory  for  the  location  of 
the  Mach  disk  is  a new  development  in  basic  theory,  and  it  is  therefore  reco-nmended 
and  requested  that  this  work  be  pushed  to  its  logical  conclusion  by  both  theoretical 
and  experimental  analyses. 
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FUTURE  WORK 


Future  work  on  this  particular  phase  of  the  blast  program  should  be  concentrated 
on  experimental  investigations  designed  to  test  the  various  hypotheses  advanced  as  to 
the  cause  of  the  Mach  disk.  Two  series  of  such  experiments  are  planned- 

A powder  charge  will  be  exploded  in  a suitably  vented  vessel  and  the  development 
in  time  of  the  jet  recorded  wi  th  high  speed  spark  or  Schlieren  photr.  graphs. 

Different  atmospheres  will  be  introduced  into  a suitably  chosen  glass  tul>e,  and  a 
jet  formed  by  exhausting  a high  pressure  (3000  psi ) tank  through  a needle  valve  whose 
opening  can  be  adiabatically  controlled,  will  be  studied. 

The  first  of  these  experiments  is  in  the  process  of  development,  while  the  second 
awaits  adequate  room  and  eouipment. 


9 


Table  I.  Correlation  of  Indices  in  Figures  I and  2 with  Numbers  of 

Characteristic  Lines 


No.  of 
Point 

C2/? 

No.  of 
Point 

V> 

C?/? 

No  .of 
Point 

V’ 

C?/? 

No. of 
Point 

c / 
/ 

1 

596 

402 

26 

594 

398 

51 

592 

394 

76 

390 

2 

597 

401 

27 

595 

397 

52 

593 

393 

77 

589.5 

390 

3 

598 

28 

596 

396 

53 

592 

393 

78 

589 

390.5 

4 

599 

399 

29 

595 

396 

r * 
Of 

591 

394 

79 

588.5 

391 

5 

598 

399 

3C 

594 

397 

55 

395 

588 

391.5 

6 

597 

31 

593 

398 

56 

589 

396 

81 

587.5 

392 

7 

596 

401 

32 

592 

399 

57 

589 

395 

82 

392.5 

8 

595 

33 

592 

398 

58 

394 

83 

586.5 

393 

9 

595 

34 

593 

397 

59 

591 

393 

84 

586.5 

392.5 

10 

596 

1 

400 

35 

694 

396 

592 

392 

85 

587 

392 

11 

597 

399 

36 

59  5 

395 

61 

591 

392 

86 

587.5 

391.5 

12 

598 

398 

37 

594 

395 

62 

393 

87 

588 

391 

13 

597 

398 

38 

593 

396 

63 

589 

394 

88 

588.5 

390.5 

14 

596 

390 

39 

592 

397 

64 

588 

395 

89 

589 

390 

15 

595 

400 

40 

591 

398 

65 

588 

394 

Q0 

} 

589.5 

389.5 

16 

594 

401 

41 

591 

397 

66 

589 

393 

91 

589 

389.5 

17 

594 

42 

592 

396 

67 

392 

92 

588.5 

390.0 

18 

595 

399 

43 

593 

395 

68 

i 

'91 

391 

93 

588 

390.5 

19 

596 

398 

44 

594 

394 

69 

590 

391 

94 

587.5 

391 

20 

597 

397 

45 

593 

394 

70 

589 

392 

95 

587 

391.5 

21 

596 

397 

46 

592 

395 

71 

588 

393 

96 

586.5 

392 

22 

595 

398 

47 

591 

396 

72 

587 

394 

97 

586.5 

391.5 

23 

594 

399 

48 

590 

397 

73 

587 

393 

98 

587 

391 

24 

593 

400 

49 

590 

396 

74 

588 

392 

99 

587.5 

390.5 

25 

593 

399 

50 

591 

395 

75 

539 

391 

100 

588 

390 
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Table  1.  Correlation  of  Indices  in  Figures  I and  2 with  Mum  ers  of 
Characteristic  Lines  (Cont'd) 


No.  of 
Point  C 


C'A 


No.  of 
Point 


Cl/l 


No.  of 
Point 


Clj'2 


No  .of 

Point  C i j-j 


534.5’ 


535 
579 

569.5 
559 
545 
533  I 

533 
535 

534 

182  534.5 

535 

535.5 
579 

569.5 
559 
545 

! .,33  ! 

i ! 

533 

I 

I 535 

I 

534  | 
534 . 5 ; 

535 
535. S1 


C ?/l 


390.5 


390.5 


390.5 
390 
390 
390 
390 
390 
390 
390 
390 
390 
390 
390 

389.5 


Table  1.  Correlation  of  Indices  in  Figures  I and  2 with  Numbers  of 
Characteristic  Lines  (Cont’d) 


Vo . o f 
Pa  i n t 

c‘/> 

Ci/ ? 

No.  of 
Point 

Ci/:> 

ct/t 

No.  of 
Point 

V’ 

C?/t 

No.  of 
Point 

V’ 

^iji 

196 

536 

389.5 

I 

221  1 

537 

388.  5 

246 

536 

387.5 

271 

533 

386.! 

197 

579 

389 

222 

579 

388 

247 

536.5 

n 

272 

533.5 

ft 

198 

569.5 

rt 

223 

569.5 

ft 

248 

537 

ft 

273 

534 

ft 

199 

559 

ft 

224 

559 

n 

249 

537.5 

ft 

274 

534.5 

ft 

200 

545 

ft 

225 

545 

n 

250 

538 

ft 

275 

535 

ft 

201 

533 

n 

226 

533 

n 

251 

579 

387 

276 

538.5 

ft 

202 

533.5 

ft 

227 

533.5 

ft 

252 

569.5 

ft 

277 

536 

»» 

2C3 

534 

- 

228 

534 

n 

253 

559 

ft 

278 

536.5 

ft 

204 

534.5 

ft 

229 

j 

534.  5 

| 

ft 

254 

545 

ft 

279 

537 

ft 

205 

535 

ft 

230 

| 535 

| 

ft 

255 

533 

ft 

280 

537.5 

H 

206 

535.5 

n 

j 231 

1 535.5 
I 

ft 

256 

533.5 

rt 

281 

538 

ft 

207 

536 

ft 

1 232 

| 

1 

536 

1 i 

ft 

257 

534 

it 

282 

538.5 

ft 

208 

536.  5 

" i 

I 233 

1 

j 536.5 

ft 

258 

534.5 

ft 

283 

539 

ft 

209 

579 

l 

i 

| 388.5  1 

234 

537 

ft 

259 

535 

ft 

284 

585 

385 

210 

569.5 

n | 

235 

537.  5 

ft 

260 

535.5 

ft 

288 

579 

•l 

211 

559 

ft 

236 

579 

387.5 

261 

j 536 

• i 

289 

509.5 

ft 

212 

545 

ft 

237 

1 

569.5 

ft 

262 

| 536.5 

ft 

290 

559 

ft 

213 

533 

' n 
1 

' 238 

559 

ft 

1 263 

| 537 

ft 

291 

545 

1 

1 W 

214 

533.5 

I " 
j 

i 239 

545 

1 " 

i 

264 

1 537.5 

1 ‘ 

292 

533 

1 " 
1 

215 

j 534 

j •» 

240 

1 533 

ft 

265 

538 

i " 

1 I 

293 

533.5 

! « 

216 

534.5 

j 

1 . 
1 

i 241 

1 533.5 

ft 

266 

538.5 

1 * 

| | 

294 

534 

217 

535 

1 n 

!!  - 

534 

ft 

j ! 

j 267 

1 

! 579 

386. 5 ' 

i 

295 

534.5 

I 

218 

535.5 

j 

i » 
1 

j|  243 

534.  5 

I " 

268 

| 569.5 

n 

| 

296 

1 

535 

385 

219 

536 

ft 

1 244 

535 

ft 

269 

559 

« j 

297 

535.5 

f» 

220 

536.  5 

1 " 

245 

535.5 

i " 

il  270 

545 

ti 

298 

536 
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Table  I.  Correlation  of  Indices  in  Figures  I and  2 with  Numbers  of 
Characteristic  Lines  (Cont'd) 


Vo  .of 
Pn  t ii  t 

C,A 

C.’  fi 

i Vo . o f 
Po  i n f 

ci/? 

o 

( *N- 

l 

No.  of 
Point 

C 1 
1/? 

C?/? 

No.  of 

Po  i n f 

V’ 

C 2/ ‘2 

299 

536.  5 

385 

326 

, 539.5 

382 

351 

559 

379 

376 

537.5 

377.5 

300 

537 

it 

327 

541 

n 

352 

545 

*t 

377 

538 

It 

301 

537.5 

n 

328 

580.  5 

380.5 

353 

533 

" 

378 

538.5 

ft 

302 

538 

n 

329 

579 

n 

354 

533.5 

" 

379 

539 

It 

303 

538.  5 

n 

330 

569.5 

n 

355 

534 

n 

380 

541 

ft 

304 

539 

n 

331 

559 

n 

356 

536 

it 

381 

544 

ft 

305 

539.  5 

n 

332 

545 

" 

357 

536.5 

n 

382 

545.5 

377.5 

305.5 

582 

382 

333 

533 

358 

537 

n 

383 

576 

376 

306 

579 

n 

o34 

533.5 

n 

359 

537.5 

it 

384 

569.5 

ft 

307 

569.5 

n 

335 

534 

n 

360 

538 

n 

385 

559 

It 

308 

559 

n 

336 

534.5 

ft 

361 

538.5 

it 

386 

545 

It 

309 

545 

n 

337 

535.5 

ft 

362 

539 

ft 

387 

533 

ft 

310 

533 

n 

338 

536 

n 

363 

541 

ft 

388 

533.5 

It 

311 

533.5 

it 

339 

536.5 

| " 

364 

544 

it 

389 

534 

ft 

312 

534 

n 

340 

537 

380.  5 

365 

1 

545.5 

| 

it 

390 

536 

It 

313 

534.5 

it 

341 

537.  5 

" 

1 

; 366 

1 

577.5 

377.5 

391 

536.5 

ft 

314 

535 

n 

342 

538 

1 

i " 

367 

j 

569.5 

It 

392 

537 

ft 

315 

535.  S 

n 

i 

343 

538.5 

1 

ft 

1 

368 

559 

It 

393 

537.5 

It 

316 

536 

| n 

344 

539 

ft 

369 

545 

It 

1 

394 

538 

ft 

317 

536.5 

it 

345 

539.5 

! " 

370 

533 

" 

| 395 

538.5 

It 

318 

537 

ft 

346 

541 

rt 

371 

533.  5 

It 

| 396 

539 

It 

319 

537.5 

r. 

347 

544 

n 

372 

i 

534 

ft 

1 397 

541 

ft 

320 

538 

n 

348 

545.5 

ri 

373 

536  i 

It 

1 

398 

544 

It 

321 

538.5 

ft 

349 

579 

379 

374 

536.5 

ft 

399 

545.5 

It 

322 

539 

n 

350  1 

1 569.5 

ft 

375 

537 

n 

403 

574.5 

374.5 
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Table  I.  Correlation  of  Indices  in  Figures  I and  2 with  Numbers  of 
Characteristic  Lines  (Cont'd) 


No  .of 
Poin . 

C?/? 

No.  of 
Point 

Ci  /? 

C?/? 

No. of 
Point 

v* 

C?/? 

No . o f 
Point 

C<A 

Cz/? 

404 

569.5 

374.5 

429 

536 

373 

454 

551.5 

371 

484 

542 

368 

405 

559 

ft 

430 

533 

ft 

455 

553 

ft 

485 

543 

ft 

406 

545 

ft 

431 

533.5 

ft 

456 

569.5 

369.5 

486 

544 

ft 

407 

533 

n 

432 

534 

ft 

457 

559 

ft 

48/ 

545.5 

ft 

408 

533.5 

ft 

433 

542 

ft 

458 

545 

ft 

488 

548.5 

ft 

409 

534 

ft 

434 

543 

ft 

459 

542 

ft 

489 

550 

ft 

410 

536 

n 

435 

544 

ft 

460 

539 

ft 

490 

551.5 

ft 

411 

536.5 

n 

436 

545.5 

f* 

461 

536 

ft 

491 

553 

ft 

412 

537 

n 

437 

548.5 

ft 

462 

541 

ft 

492 

555 

ft 

413 

537.5 

ft 

438 

550 

ft 

463 

542 

" 

493 

556.5 

ft 

414 

538 

ft 

439 

551.5 

" 

464 

| 

543 

ft 

474 

566.5 

ft 

415 

538.5 

ft 

440 

571 

371 

1 465 

544 

475 

565 

n 

416 

539 

ft 

441 

569.5 

466 

545.5 

ft 

! 476 

563.  5 

ft 

417 

541 

ft 

442 

559 

n 

i 

467 

548.5 

ft 

1 

1 

477 

562 

ft 

418 

544 

ft 

443 

545 

468 

550 

- 

| 478 

560  5 

ft 

419 

545.5 

ft 

444 

542 

n 

469 

551.5 

" 

| 494 

566.5 

366.5 

420 

542 

" 

! 

445 

539 

470 

553 

" 

495 

1 

1 

565 

ft 

421 

539 

446 

536 

, 471 

555 

ft 

496 

563.5 

h 

422 

536 

ft 

447 

541 

472 

568 

3b8 

497 

562 

ft 

423 

573 

373 

448 

542 

473 

559 

ft 

498 

560.  5 

424 

569.5 

ft 

449 

543 

479 

545 

ft 

499 

559 

425 

559 

ft 

450 

544 

480 

542 

ft 

i 

500 

545 

" 

426 

545 

r. 

451 

545.5 

„ 

481 

S39 

n 

501 

542 

ft 

427 

512 

n 

452 

548.5 

482 

536 

ft 

502 

539 

" 

428 

539 

1 « 

453 

550 

■ 

483 

541 

ft 

, 503 

536 

ft 
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Table  I.  Correlation  of  Indices  in  Figures  I and  2 with  Numbers  of 
Characteristic  Lines  (Cont'd) 


Ho. of  Ho. of  No. of  No. of 


Poin  t 

C>A 

r?/> 

Po  i n t 

r?/  7 

Point 

v* 

Cl/i 

Point 

542 

366.5 

529 

553 

365 

555 

552 

362 

550 

Em 

544 

n 

530 

555 

n 

556 

551 

581 

549 

n 

506 

545.5 

n 

531 

556.5 

o 

557 

550 

n 

582 

548 

n 

548.5 

n 

532 

558 

n 

558 

549 

ti 

583 

547 

n 

550 

n 

533 

559.5 

n 

559 

548 

n 

584 

546 

n 

551.5 

n 

535 

563.5 

363.5 

547 

n 

585 

545 

n 

553 

n 

536 

562 

n 

561 

546 

n 

586 

544 

555 

n 

537 

n 

562 

5*5 

.i 

587 

543 

n 

512 

556.5 

n 

538 

559 

n 

563 

544 

n | 

588 

542 

n 

513 

558 

n 

539 

552 

n 

564 

543 

n 

589 

541 

n 

565 

365 

540 

545 

n 

565 

542 

n 

590 

559 

359 

563.5 

n 

541 

542 

n 

566 

541 

" 

591 

558 

n 

562 

n 

542 

539 

n 

567 

n 

592 

557 

n 

560.5 

n 

543 

545.5 

n 

568 

539 

n 

593 

556 

H 

559 

H 

544 

548.5 

n 

569 

n 

594 

555 

n 

545 

n 

545 

H 

360.5 

595 

554 

n 

520 

542 

n 

546 

562 

362 

! 571 

559 

n 

596 

553 

n 

521 

539 

« 

547 

560.5 

n 

572 

558 

n 

597 

552 

n 

522 

536 

n 

548 

559 

n 

573 

557 

n 

598 

551 

n 

523 

542 

n 

549 

558 

n 

574 

556 

n 

599 

550 

n 

524 

544 

r» 

550 

557 

n 

1 

575 

555 

n 

549 

n 

525 

545.5 

ti 

551 

556 

" 

| 576 

554 

n 

548 

n 

526 

548.5 

n 

552 

555 

B 

577 

553 

n 

547 

n 

527 

550 

n 

553 

554 

B 

578 

552 

ti 

546 

n 

528 

331  . D 

» 

554 

553 

n 

579 

551 

n 

545 

tt 
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Table  1.  Correlation  of  Indices  in  Figures  I and  2 with  Numbers  of 
Characteristic  Lines  (Cont'd) 


So  .of 
Point 

r j 
1/2 

C 2/2 

No.  of 
Point 

C'  f 
Cl/2 

C2/2 

No. of 
Point 

C>/> 

C?/? 

No.  0 f 
Point 

r / 
2/2 

605 

544 

359 

630 

551 

357 

655 

553 

355 

681 

546 

356 

606 

543 

« 

631 

550 

« 

656 

552 

n 

682 

541 

359 

607 

542 

n 

632 

549 

n 

657 

551 

11 

683 

540 

360.5 

608 

558 

358 

633 

548 

n 

658 

550 

n 

684 

539 

n 

600 

557 

« 

634 

547 

» 

659 

549 

n 

685 

545 

356 

610 

556 

n 

635 

546 

n 

660 

548 

n 

686 

544 

« 

611 

555 

n 

636 

545 

tt 

661 

547 

rt 

687A 

543 

ft 

612 

554 

n 

637 

544 

ft 

662 

554 

354 

687 

542 

358 

613 

553 

ft 

638 

543 

»i 

663 

553 

ft 

688 

546 

355 

614 

552 

ft 

639 

n 

664 

552 

n 

689 

545 

fl 

615 

551 

« 

640 

ft 

665 

551 

n 

690 

552 

352 

616 

550 

ft 

ft 

666 

550 

ft 

691 

551 

ft 

617 

549 

ft 

ft 

667 

549 

ft 

692 

550 

ft 

618 

548 

ft 

556 

356 

668 

548 

n 

693 

549 

ft 

619 

547 

r 

555 

ft 

669 

547 

ft 

694 

548 

ft 

620 

546 

ft 

645 

554 

n 

670 

546 

ft 

695 

547 

n 

621 

545 

n 

646 

553 

n 

671 

545 

ft 

696 

551 

351 

622 

544 

f» 

647 

552 

n 

672 

553 

353 

697 

550 

« 

623 

543 

ft 

648 

551 

n 

673 

552 

ft 

698 

549 

624 

557 

357 

649 

550 

n 

674 

551 

tt 

699 

1 

550 

350 

625 

556 

f» 

650 

549 

n 

675 

550 

n 

700 

549 

« 

626 

555 

n 

651 

548 

n 

676 

549 

tt 

703 

551.5 

363.  5 

627 

554 

n 

652 

54/ 

n 

677 

548 

n 

704 

553 

tl 

628 

553 

i» 

653 

555 

355 

678 

547 

n 

505 

555 

n 

629 

552 

ti 

654 

554 

n 

679 

706  | 

556.5 

n 

16 


Table  I 


. Correlation  of  Indices  in  Figures  I and  2 with  Numbers  of 
Characteristic  Lines  (Cont'd) 


No.  of 
Point 

c</» 

C?/? 

No. of 
Point 

W> 

*"?/ ? 

No. of 
Point 

c‘/> 

Cl/l 

No. of 
Poin  t 

V’ 

C»/« 

558 

363.5 

732 

558 

359 

747 

568 

357 

772 

565.5 

354 

559.5 

n 

733 

559.5 

ft 

748 

569 

It 

773 

567 

ft 

709 

561 

ft 

734 

561 

ft 

749 

556.5 

356 

774 

568 

It 

551.5 

362 

735 

562.5 

ft 

558 

n 

775 

569 

ft 

711 

553 

n 

736 

564 

ft 

751 

559.5 

ft 

776 

570 

ft 

712 

555 

it 

737 

565.5 

It 

752 

561 

ft 

777 

571 

ft 

713 

556.5 

it 

738 

567 

It 

753 

562.5 

it 

678 

547 

353 

714 

558 

n 

688A 

556.5 

358 

754 

564 

■ 

778 

546 

It 

715 

559.5 

n 

689A 

558 

n 

755 

565.5 

it 

779 

545 

It 

716 

561 

ft 

690A 

559.5 

it 

756 

567 

n 

780 

561 

ft 

717 

it 

691 A 

561 

it 

757 

568 

it 

781 

562.5 

It 

718 

551.5 

360.5 

692A 

562.5 

ft 

758 

569 

it 

782 

565.5 

ft 

719 

553 

it 

693A 

564 

ft 

759 

544 

355 

783 

567 

ft 

720 

555 

n 

694A 

565.5 

ft 

760 

559.5 

ft 

784 

568 

It 

721 

556.5 

• 

695A 

567 

ft 

761 

561 

It 

785 

546 

352 

722 

558 

n 

696  A 

568 

ft 

762 

562.5 

It 

786 

545 

ft 

723 

559.5 

n 

639 

542 

357 

763 

564 

n 

787 

548 

361 

724 

561 

*t 

739 

556.5 

ft 

764 

565.5 

« 

1 

7(,8 

547 

ft 

725 

562.5 

n 

740 

558 

n 

765 

567 

h 

789 

546 

ft 

545 

l» 

726 

564 

ft 

741 

559.5 

ft 

766 

568 

f» 

790 

561 

ft 

727 

565.5 

ft 

742 

561 

It 

767 

569 

ft 

791 

549 

350 

728 

551.5 

359 

743 

562.  5 

ft 

768 

545 

354 

I'iZ 

548 

ft 

729 

553 

B 

744 

564 

" 

769 

561 

ft 

793 

547 

r> 

730 

555 

ft 

745 

565.5 

It 

770 

562.5 

ft 

794 

546 

ft 

731 

556.5 

ft 

746 

567 

ft 

771 

564 

ft 

795 

545 

ft 
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Table  1.  Correlation  of  Indices  in  Figures  I and  2 with  Numbers  of 
Characteristic  Lines  (Cont'd) 


No. of 
Point 

V* 

No. of 
Point 

f'  / 

' 1/2 

^2/2 

No. of 
Point 

V; 

<y> 

No.  of 
Point 

Cl/ 2 

796 

544 

350 

807 

562.5 

352 

818 

568 

351 

828 

569 

797 

543 

ft 

808 

565.5 

ft 

819 

569 

ft 

829 

798 

542 

ft 

809 

567 

ft 

820 

570 

ft 

547 

561 

799 

549 

349 

810 

568 

ft 

821 

571 

« 

831 

565.5 

800 

548 

ft 

811 

569 

ft 

822 

572 

ft 

832 

567 

801 

570 

355 

812 

570 

ft 

823 

561 

833 

568 

802 

571 

ft 

813 

571 

ft 

824 

562.5 

ft 

834 

569 

803 

569 

353 

814 

572 

It 

825 

565.5 

ft 

835 

804 

570 

ft 

815 

562.5 

351 

826 

567 

ft 

836 

548 

805 

571 

ft 

816 

565.5 

ft 

827 

568 

ft 

837 

567 

806 

572 

ft 

817 

567 

ft 

I 


i 

i 

I 
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r2  /? 
350 

ft 

349 


348 

365.5 


Table  II.  Numerical  Relations  for  the  Two-d imons ional  Characteristics  System 

The  table  Rives  as  a function  of  the  characteristic  number  N (column  1)  the  following 
auant i ties: 

Column  2 The  value  of  v,  deviation  of  expansion  from  Mach  numlier  1. 

Column  3 The  value  of  the  stream  Mach  number  (determined  from  column  2). 

Column  4 The  ratio  of  the  local  static  pressure  p to  the  critical  static  pressure  p* 


corresponding  to  the  pressure  for 
Column  5 The  Mach  angle  u.  determined  from 
Note:  ^ l + ^2  = N 


2 

/»= 


1 + C2 
2 

V 

H 

P 

P* 

u. 

1000 

0° 

1.0000 

1 . 00000 

90° 

999 

i° 

1.0813 

.90724 

67°37 ' 

998 

2° 

1.1327 

.85100 

61=59' 

997 

3° 

1.1764 

.80486 

58° 13  ‘ 

996 

4 = 

1.2185 

.76192 

55°  9' 

995 

5° 

1.2577 

.72329 

52=39 ' 

994 

6° 

1. 2938 

. 68896 

50°37 ' 

993 

O 

/ 

1 . 3298 

.65588 

48°45' 

992 

8° 

1 . 3648 

. 62488 

47=  6' 

991 

9° 

1 . 4000 

.59484 

45=35' 

990 

10° 

1.4350 

. 56608 

44" 10' 

989 

11° 

1 . 4690 

.53921 

42°54 ' 

988 

12° 

1.5032 

.51325 

41°42 ' 

987 

13° 

1 . .368 

. 48875 

40c36 " 

986 

14° 

1.5706 

.46512 

39° 3.'  ' 

985 

15° 

1 . 6046 

.44233 

38°33 ' 

984 

16° 

1.6382 

.42075 

37=37 ' 

983 

17° 

1.6721 

. 39995 

36 "44' 

(M  = 1). 
olumn  3. 


H 

V 

P 

P* 

- 

963 

37° 

2.4107 

. 12735 

24=30' 

962 

38 

2.4525 

.11970 

24=  4' 

961 

39" 

2.4946 

11177 

23=37' 

960 

40” 

2.5373 

. 10454 

23=12' 

959 

41” 

2.5812 

.09770 

22=48' 

958 

42° 

2.6256 

.09117 

22=23' 

957 

43" 

2.6712 

. 08500 

21=59' 

956 

44“ 

2.7175 

.07917 

21=36' 

955 

45° 

2.7650 

.07359 

21=13' 

954 

46° 

2.8118 

. 06848 

20'50' 

953 

47° 

2.8601 

. 06363 

20=28' 

952 

48° 

2.9112 

. 05895 

20°  5' 

951 

49  = 

2.9613 

. 05466 

19=44' 

950 

50= 

3.0122 

. 05065 

19=23' 

949 

51  = 

3.0650 

. 04660 

19=  2' 

948 

52  = 

3.1188 

.04310 

18=42' 

947 

53° 

3.1735 

. 03985 

18=22' 

946 

54° 

3.2263 

. 03669 

18°  3' 
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Table  II.  Numerical  Relations  for  the  Two-dimensional 
Characteristics  System  (Cont'd) 


Nr. 

'l  *C2 
2 

V 

V 

P 

P* 

u 

V 

V 

¥ 

P 

P* 

a 

982 

18° 

1. 7064 

.37981 

35°53' 

945 

55° 

3.2866 

.03375 

17°43 ' 

981 

19° 

1. 7408 

.36053 

35°  4' 

944 

56’ 

3.3455 

.03103 

17°24 ' 

980 

20° 

1.7753 

.34210 

34° 16' 

943 

57° 

3.4065 

. 02843 

17°  4' 

979 

21° 

1.8098 

.32454 

33°32 ' 

942 

58" 

3. 4682 

. UZ0U5 

16"45 ’ 

978 

22° 

1.8446 

. 30768 

32°49 ' 

941 

59° 

3.5300 

.02382 

16°27 ' 

977 

23° 

1.8796 

.29152 

32°  8' 

940 

60° 

3.5945 

.02175 

16°  9' 

976 

24° 

1.9148 

. 27610 

31 °29 ' 

939 

61° 

3.6605 

.01984 

i5°sr 

975 

25° 

1.9503 

.26136 

30°50 ' 

938 

62° 

3.7283 

.01838 

15°33 ' 

974 

26° 

1 . 9860 

.24727 

30° 14 ' 

937 

63° 

3.7975 

.01643 

15°16' 

973 

27° 

2.0235 

. 233  24 

29°37 ' 

936 

64° 

3.8699 

.01489 

14°58 ' 

972 

28° 

2.0600 

. 22034 

29°  2' 

935 

65° 

3.9425 

.01350 

14°42 ' 

971 

29° 

2.0968 

. 20805 

28°29 ' 

934 

66° 

4.0175 

.01223 

14°25 ' 

970 

30° 

2. 1320 

. 19690 

27°58 ' 

933 

67° 

4.0950 

.01103 

14°  8' 

969 

31° 

2. 1726 

. 18483 

27°24 ' 

932 

68° 

4.175 

. 00990 

13°56 ' 

968 

32° 

2.2105 

. 17430 

26°54 ' 

931 

i 

69° 

4.256 

.00889 

13°35 ' 

967 

33° 

2. 2485 

. 16400 

26°24 ' 

930 

70° 

4.  340 

.00800 

13° 19 ' 

966 

34° 

2.2879 

. 15406 

25c55 ' 

929 

71° 

4.427 

.00720 

13°  3' 

965 

35° 

2.3290 

. 14480 

25°25 ' 

928 

72° 

4.517 

. 00646 

12°47 ' 

964 

36" 

2.3693 

. 13582 

24°58 ' 

927 

73° 

4.610 

.00577 

12°32 ' 

926 

74° 

4.704 

.00510 

12° 16 ' 

909 

91° 

7.008 

. 000455 

8°  7' 

925 

75° 

4.802 

. 00450 

12°  r 

908 

92° 

7.202 

.000381 

7°  54' 

924 

76 c 

4.905 

. 00400 

11°46' 

907 

93° 

7.407 

.000318 

7"41 ' 

923 

77° 

5.O10 

. 00354 

1 1°31 ' 

906 

94° 

7.622 

.000265 

7°  28' 

Q?2 

78° 

5.  119 

. 003 13 

1 1 0 16 ' 

905 

95° 

7.853 

. 000220 

r is' 
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Table  II.  Numerical  Relations  for  the  Two-d.mens ional 
Characteristics  System  (Cont’d) 


/V= 


c,  + s 

V 

M 

n 

o* 

U 

V 

V 

* 

P 

P * 

u 

2 

921 

79° 

5.231 

.00277 

n°  r 

904 

96” 

8.092 

.000182 

7°  3' 

920 

80° 

5.348 

.002413 

10c47' 

903 

97° 

8.343 

.000149 

6°50' 

919 

81 c 

5.472 

.002106 

10°32 ' 

902 

98° 

8.616 

.0001194 

6°37 ' 

918 

82° 

5.599 

.001835 

10° 17  ' 

901 

99* 

8.902 

. 0000965 

6°25' 

917 

83° 

5.  730 

.001598 

10°  3' 

900 

100n 

9.210 

. 0000770 

6°12 ' 

916 

84° 

5.865 

.001383 

9°41 ' 

899 

101  c 

9.535 

.0000614 

5 ”59 ' 

915 

O 

oc 

6.C07 

.001193 

9°27' 

898 

102'5 

9.881 

.0000481 

5°47  * 

914 

86° 

6.  152 

.001023 

9°14' 

897 

103° 

10.260 

.0000377 

5°34 ' 

913 

87° 

6.  307 

.000879 

9"  1 ' 

896 

104° 

10.665 

. 0000294 

5n21 ' 

912 

88° 

6.472 

.000751 

8°47' 

895 

105° 

11.088 

.0000224 

5°  9' 

911 

89° 

6.642 

. 000638 

8°34' 

894 

106c 

11.552 

. 00C0168 

4°57' 

9i0 

90° 

6.821 

. 000540 

8=20'  | 
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APPENDIX  A 


The  Method  of  Characteristics  for  two-dimensional,  steady,  isentropie  flow  in  two 
dimensi ons . 


I . Theory 

Let  u and  v be  the  x and  y components  of  velocity  q,  respectively,  and  let  U oe 
the  angle  the  vector  velocity  q makes  with  some  fixed  reference  line  in  the  x,  y plane. 
Let  p he  the  pressure,  C specific  heat  at  constant  pressure,  C - specific  heat  at 

V c’ 

fi 

constant  volume,  r the  density,  k - the  ratio  of  specific  heats  q~.  The  differential 

equations  of  motion  for  steady,  perfect,  i r rot  at iona 1 , isentropic  flow  in  two  dimen- 
sions are:  (see,  e.  g.  Courant  and  Friedrichs,  Ix>c . Cit.) 


-1 

7(P*) 


uv  ♦ vv 


-1 

7tPv> 


(1) 


(ru),  * (rv)y  - 0 


(2) 


v - u - 0 (3) 

x y v ' 

p - p(  r ) only 
p ■ ark 


a cons  t . 


From  the  sound  speed  c,  c* 
(2)  gives 


iiE,  we  have  r -*  — px.  Tliis  together  with  (1)  and 

r C 


(u2  - c^u,  ♦ uv(uy  * vj  ♦ (v2  - c2)vy  - 0 


(4) 


The  relation  (3)  implies  the  existence  of  a function  q?  (x,y),  called  the  potential 
function,  which  satisfies 


u * vx,  v = ^y 

In  (4)  this  gives: 

(u2  - c2)  oxx  * 2u vQpxy*(v2  - c2)9nv  = 0.  (5) 

a nonlinear  differential  equation  for  the  determination  of  the  potential  9 (x.y).  It 
is  wtll  known  that  such  a partial  differential  equation  defines  a family  of  curves 
S(x.y)  ■=  Const.,  called  the  characteristics.  They  are  obtained  by  the  substitution 
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S 2 = ? 


Then  from  (5): 


, S 2 = T , S S - <p  . 
y yy*  * y r*y 

(u2  - c2)  Sx2  + 2uvSxSr  + (v2  - c2)Sy2  = 0 (6) 

dy 

The  slope  of  a surface  S(x,y)  is  given  by  m = — - — so  that  (6)  gives  (u2  - v2)m2  - 

dx  S 

y 

2uvm  (v2  - c2)  = 0.  This  quadratic  equation  has  the  two  solutions 


m * | 
m-  \ " 


2uv  +_\j 


4ii2v2  - 4(u2  - 


2)  (v2 


C2) 


2(u2  - v2) 


Evidently,  at  any  point  in  the  physical  (x,y)  plane  there  are  defined  two  curves,  de- 
pending on  velocity  solutions  u,v  of  equation  (5).  These  slopes  are  real  if  and  only 
if 

o n n I)  I) 

4u^vz  - 4(u*  - c*  ) (v*-  - c*-)  >0 

i.e.,u2-»v2=q2^.c2  (7) 

the  flow  is  supersonic.  If  the  above  quantity  (7)  is  equal  to  0 at  a point,  the  flow 
is  sonic  there,  and  there  is  only  one  characteristic  at  the  point.  Evidently  the  net 
formed  by  the  two  characteristics  through  each  point  is  degenerate  on  r.  sonic  line  in 
the  flow. 

The  characteristics  can  lie  shown  to  have  the  following  properties: 

1.  Discontinuities  in  the  derivatives  of  u and  v can  only  occur  on  these 
characteristic,  or  "Mach",  lines.  Thus  the  Mach  lines  can  be  visualized  as  wave  fronts 
in  supersonic  flow. 

2.  The  characteristics  at  a point  make  equal  angles  A with  the  streamline, 
and  this  angle  is  the  so-called  Mach  angle.  We  find 

. c 
sin  A - — 

q 

3.  The  characteristic  net  represents  a coordinate  system  in  the  physical 
plane  (except  on  degenerate  sonic  lines).  If  we  introduce  parameters  a,  b along  the 
Mach  lines,  we  have  the  characteristics  given  by 

C + : = m.x. 

c* : yb  = m-*b 
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We  may  interchange  dependent  and  independent  variables  in  enu.it  ion  (4)  by 
the  transformation 


",  = jy,-  l«y  )*,•  v,  - ivu*  vy  - ix,. 

as  long  as  the  Jacobian  determinant  j - u>Vy  - uyvx  tloes  not  vanish  in  the  r(.gion  ^der 
consideration.  We  find,  analogous  to  eouations  (3)  and  (4), 

(a)  - yu  - 0 

(b)  (u2  - c 2 )y  y - uv(xt  + yu)  + (v2  - C2 )xu  = 0 (81 

By  the  first,  there  exists  a function  (u,v),  a "hodograph  potential,"  which  satisfies 

(u2  - c2)  4>yy  - 2uv  4>ux  * (v2  - c2)  4>uu  - 0 (9) 

v (u,v)  is  the  Legendre  (contact)  transformation  of  !P(x,y):  f>  - ux  * vy  - 3>.  Note  that 

this  eauation  is  linear.  A new  set  of  characteristics  T(u,v)  is  found  as  before  to 
satisfy: 


(u2  - c2)Ty2  - 2uvT(,Ty  i (v2  - c2)Tu2  = 0 


dll  ~ ^ 

The  slope  of  these  u,  v characteristics  is  found  to  lie  — = — aact 

dv  -m  + 


G * : ti  = - m v 

» • • 

G*:  ufa  - -m,vfc  (10) 

1.  These  characteristics  are  independent  of  the  solution  x - x(u,v.),  y - y(u,v), 
and  are  thus  fixed  in  the  u,v  plane. 

2.  The  relations  uyxb  * vxyfc  = 0 can  be  obtained,  and  indicate  the  following 

"bx.  + vb>’.  = 0 

vital  fact:  Consider  the  image  point  u,v  of  a physical  point  x,y  (j(x.y)  = 0).  Tlie 

C*  characteristic  through  x,y  is  perpendicular  to  the  C"  characteristic  through  u.v, 
and  the  C*  characteristic  through  x,.y  is  perpendicular  to  the  G*  character!  st  i c through 
u.v.  This  is  the  basis  of  the  graphical  construction  method.  Eouations  (10)  are  in- 
tegrated directly  using  Bernoulli’s  equation  for  c: 


u 


V 


(1  - u. 2 )c 


, 2 


k - 1 
k ♦ 1 
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One  finds  the  characteristics: 


\ 

/ k + 1 

q2 

* - N 

* Sin ‘ 1 

k - 1 

(k-1)  — - k 

_ <~*2  _ 

(It) 


♦ Sin ' 1 


(k 


* 2c 


wfiere  q = c*  at  M = 1.  These  curves  are  epicycloids,  and  are  generated  by  the  motion 

of  a rim  point  of  a circle  of  diameter  c*  (—  - 1)  rolling  on  a circle  of  radius  c*  . 

u 

The  pressure  ratio  is  related  to  the  velocity  by 


q 2 k-1  k*l 

— 2 , (1  - p/Po)-  (—  ) (12) 


where  po  is  the  stagnation  pressure  (chamber  pressure).  Thus  the  pressure  is  constant 
on  circles  concentric  with  the  sonic  circle  in  the  u,v  plane. 

For  practical  use,  it  is  better  to  consider  ')  positive.  This  can  be  done  by  tak- 
ing '3  = 200  degrees  for  an  undisturbed  flow,  v = 0.  (See  p 36,  A.  Ferri.)  Every 
characteristic  can  then  be  distinguished  by  the  value  of  the  constant  C in  equation 
(11).  For  the  G+  characteristic. 


-0  - f (q  ) - c , (13) 

where  C(  is  the  absolute  magnitude  of  the  constant  that  in  equation  (11)  is  negative. 
For  the  G~  characteristic, 


<:  = f(q)  - c 2 ; 

eliminating  f,  then  0 in  (i3)  and  (14),  we  have 


0 = 


f(q) 


(14) 


(15) 
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From  (12)  and  (15),  f(o ) is  determined  by 


2 P 

— and  the  cor  responding  values  of  — , 

2 p* 

the  Mach  number  M,  the  Mach  angle,  and  the  angle  for  the  corresponding  Prandtl -Meyer 


flow  are  given  in  Table  II.  Thus,  entering  in  Table  VI  the  values  of 


2 


for 


points  in  the  f'ow,  taken  from  Table  I,  yields  the  values  of  the  flow  quantities  at 
tliose  points. 


II.  Application  of  method  of  characteristics  to  the  construction  of  the  flow  through 

a nozzle. 

Given  anozzle  symmetric  about  the  center  line  (which  is  chosen  to  lx*  the  x-axis), 
locate  the  throat  of  the  nozzle  at  the  y-axis.  Now  plot  the  nozzle  contour  so  that  it 
is  congruent  vith  the  actual  contour  but  of  convenient  size.  The  hodograph  is  now 
mounted  so  that  the  line  whose  direction  numlier  is  200  is  paral led  to  the  x-axis. 

Since  the  Mach  net  is  degenerate  on  the  sonic  line  x - 0,  choose  the  initial 

curve  in  this  way:  Through  that  point  where  the  nozzle  has  at  least  a slope  of  1°  (or 

a multiple  thereof)  pass  a circular  arc  wliose  center  is  the  intersection  of  the  tan- 
gent to  the  nozzle  contour  at  that  point  with  the  x-axis.  By  the  hydraulic  theory, 
compute  the  flow  Quantities  on  this  arc.  From  Table  II,  find  the  pressure  nimber,  N. 
Divide  the  arc  into  equal  sub-arcs  whose  included  angle  is  1’  (or  a multiple  thereof). 
At  the  end  point  of  each  arc  label  the  point  in  a suitable  fashion.  The  pressure 
number  corresponding  to  each  point  is  N.  The  direction  number  is  200  plus  the  in- 
cluded angle  (in  degrees)  with  the  x-axis.  The  two  characteristics  in  the  hodograph 
plane  which  pass  through  the  image  of  each  point  can  now  be  deduced.  Two  adjacent 
points  in  the  x,y  pl^ne  determine  a third  point  as  follows: 

Let  the  upper  point  be  A,  the  lower  one,  B.  Corresponding  to  the  upward-going  Gf 
through  the  image  of  A in  the  hodograph  plane,  there  is  a downward  sloping  Mach  line 
(G' ) through  A,  perpendicular  to  the  G*  characteristic.  Similarly,  there  is  a line  C, 
through  B,  sloping  upwards,  perpendicular  to  the  G"  (downward  sloping)  characteristic 
through  the  image  of  B in  the  hodograph  plane.  The  intersection  of  the  C"  line  from 
A and  the  CA  line  from.  B determines  a third  point,  P,  where  the  flow  Quantities  are 
determined  by  the  location  of  the  intersection  of  the  co r respond ing  G*  and  G"  charac- 
teristics in  the  hodograph  plane.  This  process  may  now  be  continued. 

At  boundaries,  if  the  boundary  is  a rigid  wall,  the  new  point  is  determined  by 
one  characteristic  and  the  fact  that  the  direction  numlier  must  correspond  to  the  slope 
of  the  wall.  At  a free  surface,  the  new  point  is  determined  by  one  characteristic  and 
the  requirement  that  the  pressure  index,  N,  has  a certain  value  which  corresponds  to 
receiver  pressure. 
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III.  Construction  of  the  Ho'ioiiraph 


Tin'  equipment  usually  used  includes  a hodograph  with  a suitable  ellipse  and  a 
drafting  tool.  Their  use  is  outlined  in  I.  From  the  knowledge  of  the  pressure  and 
direction  of  flow  at  a point,  one  can,  using  the  hodograph.  determine  the  Mach  angle 
and  the  di-ections  of  the  characteristics  passing  through,  that  point.  Or,  knowing 
which  two  characteristics  pass  through  a point  in  space  one  can  find  the  flow  Quanti- 
ties there  by  using  the  hodograph. 

In  order  to  fix  the  quantities  involved  in  such  calculations,  it  is  necessary  to 
index  the  characteristics  in  a standard  fashion.  That  method  of  indexing  used  here  is 
the  one  proposed  by  Buseman  and  is  in  general  use.  The  procedure  is  as  follows: 

1.  Choose  a radius  for  the  sonic  circle. 

2.  The  Ixidograph  can  now  be  drawn  according  to  the  directions  given  in  I. 

3.  Draw  one  characteristic  leaving  the  sonic  circle  every  two  degrees,  curving 
downward,  and  another  set,  leaving  from  the  same  points,  curving  upward. 

4.  Number  the  lines  in  each  set  so  that  at  any  point  along  the  sonic  circle  the 
sum  fa)  of  the  number  given  one  line,  and  (b)  of  the  number  f run,  the  line  at  that 
point  which  he  longs  to  the  other  family  is  1000. 

5.  At  that  point  which  is  the  in te r sec f ion  of  the  sonic  circle  with  the  radius 
which  is  paralleled  to  the  arbitrarily  chosen  fixed  line  in  space,  make  the  difference 
of  the  numbers  referred  to  in  (4)  be  200.  Then  the  sum  (called  the  pressure  index)  of 
the  numljers  of  the  two  characteristics  through  some  point  in  space  gives  the  ratio: 
Critical  speed/flow  speed.  (This  relationship  is  shown  in  Table  II.)  The  difference 
of  the  two  numbers  gives  the  direction  of  the  flow  with  respect  to  the  axis  fixed  in 
space.  From  the  former  number  other  pertinent  quantities  can  be  obtained  (Table  II). 

The  initial  curve  is  taken  to  be  the  arc  of  the  circle  whose  center  is  the  vertex 
of  the  norzle  and  which  passes  through  the  endpoints  of  the  throat  section.  On  that 
line  it  is  assumed  that  the  flow  ve  loci ty- loca 1 sound  velocity  (eounls  throat,  or 
critical,  velocity)  and  the  pressure  index  is  then  1000  there.  The  flow  is  now  de- 
termined by  fixing  the  atmospheric  pressure,  which  is  assumed  to  be  14.7  ps i . 
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Figure  1.  Flow  through  a nozzle  of  7°  half-angle 
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Figure  3.  Pressure  along  the  axis  vs  distance  from  throat 
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Figure  4.  A typical  shock  bottle  formed  at  a gun  muzzle 
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